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ABSTRACT 

There are by now ten published detections of fast radio bursts (FRBs)—single bright 
GHz-band millisecond pulses of unknown origin. Proposed explanations cover a broad 
range from exotic processes at cosmological distances to atmospheric and terrestrial 
sources. Loeb, Maoz, and Shvartzvald have previously suggested that FRB sources 
could be nearby flare stars, and pointed out the presence of a W-UMa-type contact 
binary within the beam of one out of three FRB fields that they examined. To fur¬ 
ther test the flare-star hypothesis, we use time-domain optical photometry and spec¬ 
troscopy, and now hnd possible flare stars in additional FRB fields, with one to three 
such cases among all eight FRB fields studied. We evaluate the chance probabilities of 
these possible associations to be in the range ^ 0.1% to 9%, depending on the input 
assumptions. Further, we re-analyze the probability that two FRBs recently discovered 
3 years apart within the same radio beam are unrelated. Contrary to other claims, we 
conclude with 99% confidence that the two events are from the same repeating source. 
The different dispersion measures between the two bursts then rule out a cosmological 
intergalactic-medium origin for the dispersion measure, but are consistent with the 
flare-star scenario with a varying plasma blanket between bursts. Finally, we review 
some theoretical objections that have been raised against a local flare-star FRB origin, 
and show that they are incorrect. 
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1 INTRODUCTION 


Nearly a decade after their discovery by iLorimer et al.l 
(120011 ). the source and nature of fast radio bursts (FRBs) 
is yet unknown. FRBs are bright 0.1 — 1 Jy) and brief 
('^ 1 ms) pulses of ~ 1 GHz radio emission. The latest es¬ 
timate of the rate of FRBs with flux > 0. 1 Jy is 3.31^ x 


ten candidate bursts have been reported ( 

--, 

Lorimer et al. 

2007 

Keane et al. 20121: Thornton et all 201f 

: Snitler et al. 

2014 

Burke-Snolaor & Bannisteill20l4lPetroff et al.ll2015al: 

Ravi et all 2015|l. Nine of them have been discovered with 


the Parkes Radio Telescope, and one with the Arecibo Radio 


Telescope. Eight of the FRBs were found in the analysis of 
archival data, long after the e vents. Two FRBs were discov¬ 
ered in real time (within 10 s) llRavi et aI.l2015l : IPetroff et al.l 
l2015al) . one of them while observing the Carina dwarf galaxy, 
the other while re -monitoring the locati on of a previous 
burst (FRB110220. [riiorntonA all (l2013l )). The latter was 
also the first FRB with a polarization measurement, and was 
found to be up to ~ 50% circularly polarised, depending on 
the temporal part of the pulse considered. In no cases have 
there been known transient or eruptive counterparts at other 
wavelengths coinci dent with the FRB, i ncluding in the real¬ 
time detection bv IPetroff et all ll2015al ). where prompt fol¬ 
lowup observations at other bands were undertaken (within 


© 0000 RAS 







































2 D. Maoz et al. 


9 hrs in X-rays and within 17 hrs in the optical and near 
infrared). No clear consensus has emerged on whether FRBs 
are Isotropically distributed, or avoid the Galactic plane, due 
to some combination of ba ckgrounds, selection effects, and 
prop a gation effects (see, e.g. Petroff et al ] |2015al : [Rane et all 


1201^ : iMacauart fc Johnstori l2015l l. An intrinsic avoidance 
of the Galactic plane is not expected in any model that has 
been proposed for FRBs. 

The dispersion measures (DMs) of FRBs, i.e. the line- 
of-sight column densities of free electrons, as deduced from 
the normalisation of the f~^ frequency dependence of the 
pulse arrival times, exceed the values expected from models 
of the Milky Way interstellar electron distribution. FRBs 
have DMs in the range ~ 300 — 1600 pc cm~^, or elec¬ 
tron column densities Ne ~ (1 — 5) x 10^^ cm“^. FRBs 
have therefore been inferred to originate from extragalactic 
sources at cosmological distances. The intergalactic medium, 
perhaps combined with the DM of the sources’ host galax¬ 
ies, would account fo r the excess DM va lues. In one case, 
however fFRB010621. lKeane et ahl (l2012l ~)l. direct estimates 
of the ionized gas column density in this direction show 
that the Galactic models underestimate the DM, placing 
the FRB source within the Galaxy at a distance of 8-20 kpc 
dBannister fc Madsen|[20l3 l, or possibly much closer If some 
of the DM is intrinsic to the source. In addition to the ex - 
cess DM , three or four o f the F RBs (iL orirner et al.1 ll2007ll . 
one from Thorntoii_et_aljj201^ Ravj_et al.1 ( 2015ll . and pos¬ 
sibly iBurke-Spolaor fc Bannister ( 2014l ll show frequency- 
dependent pulse smearing, indicative of scattering by an in¬ 
homogeneous electron distribution along the line of sight, at 
a level much larger than expected by the Galactic interstellar 
medium. The intergalactic medi um has been invoked to e x¬ 
plai n this scattering a s well fe.g. IMacauart fc Koavl ll2013ll l. 
but iKatj (I2OI4I . [2OI5I I has argued that the implied densities 
are too high and the scattering instead likely occurs near 
the source. 

If FRB sources are extragalactic, their bright fluxes 
imply high isotropic luminosities of Liso ^ 10"^^ erg s“^. 
A flurry of theoretical ideas has been put fo rward to ex¬ 
plain them, including: black- hole evaporation (iKeane et al. 
2012); magnetar hyperflares (|Po£0^.^££ostnOT 201^ 


201T); neutron-star mergers (iTotanif 2013 Ravi fc Laskv 
2014); white-dwarf mergers llKashivama et ah 20]j); c^ 


lapse of supram assive neutron stars ( Falcke fc Rezzolla 
[201^ Zhanj2^M ); orbi ting bodies immersed in pulsar winds 
( Mottez_&Zaxlra [203)j magnetar pulse-wind interactions 
( Lyubarsky 2014 ); giant pulses from magnetars near galac¬ 
tic centers ( Pen fc C onnoill2 015ll: collision s between neutron 
stars and asteroids (iGeng fc Huandl2015l '): giant pulses from 
you ng pulsars within su pernova remnants in nearby galax¬ 
ies llConnor et aDl2015h or at cosmolog ical distances ( Kat j 
l2015[ l: quark novae l ^and et al.|[2015l l ; and even directed 
signa ls from extraterrestrial civilizations llLuan fc Goldreid] 
l2014h . 

In parallel, the possibility has lingered that 
FRBs have a terrestrial, possibly human-m ade, origin 
llBurke-Spolaor et al.1120111 : iKulkarni et al]|2014l l. From the 
start, some resemblances were noted between FRBs and 
“perytons” llBurke-Spolaor et al.ll20iTl l. signals with fluxes, 
time structures, and roughly quadratic frequency sweeps 
similar to FRBs but, as opposed to FRBs, with detections 
in all focal-plane receivers of the telecope, indicating a 


nearb y source that has not been focused by the antenna 
dish. IPetroff et al.1 ll2015bl l have recently demonstrated 
that perytons arise from microwave kitchen ovens at the 
Parkes observatory (generally during lunchtime), when an 
oven door is opened while the appliance is in operation. 
They presented further statistical evidence that the FRBs 
are a truly astronomical phenomenon, distinct from the 
human-made perytons. 

Alternatively to thes e exp lanations for FRBs, 
iLoeb. Shvartzvaid. fc Maoz I ||2014 1 proposed that FRB 
sources could be nearby flaring stars. Abundant types 
of flare stars that might provide the required FRB rates 
include young M-type stars, W-UMa-class stars (stable 
contact binaries composed of two stars of about solar 
mass), or other types. Flaring main-sequence stars are 
already known to produce coherent radio bursts with 
millisecond-scale d urations and ~ 1 GHz flux densities of 
a fraction of a Jy (iLang et al.lll983l: iLang fc Willsonlll986l : 
iGudel et al.1 1 19891 : Bastian et al. ~ ll990l l. Although details 
of the flare emission mechanisms are poorly understood, 
these bursts are thought to be produced by a cyclotron 
maser process, in which bunches of electrons emit coherent 
radio emission as they gyrate around the ma g netic field 
of the host star (iGiidell I2OO2I : FiV euman 3 12OO6I : lMatthev3 
I2OI3II . iLoeb et al. I 2014h argued that if FRB eruptions 


are produced at the bottom of the coronae of their host 
stars, or the erupt i ons c ause coronal mass ejections of their 
own JPrake et al.1 l2013l l. the radiation will pass through 
a plasma blanket with a characteristic electron column 
density > 10^° cm“® x Rq ~ 300 cm“® pc, as needed to 
explain the observed DMs of FRBs. Such co lumn densities 
are typ ically inferred for stellar flares (iGetman et al.1 
l2008allal We note that radio bursts from stellar flares are 
observed to be circularly polariz ed, like the single FRB 
with a polarization measurement (jPetroff^^^l^[20153). 

T o test the flare-star FRB hypothesis, iLoeb et al. I 
(I2OI4II monitored in the optical the fields of three FRBs for 
variable stars. In one field, they discovered a nearby (800 pc) 
W UMa system at a position within th e FRB radio beam . 
From the abundance of W UMa systems. iLoeb et al. I (l2014h 
derived a probability of ~ 5% for a chance coincidence of 
such a system in the beam of one out of three FRB beams. 
However, this inference was compromised by the a posteri¬ 
ori nature of the conclusion, since the objective of specifi¬ 
cally finding W UMa stars was not defined in advance of the 
search. 

In this paper, we extend the observational test of the 
flare-star hypothesis, by means of photometric and spec¬ 
troscopic monitoring and analysis of stars in 8 out of 9 
known FRB fields. We find one new example of a W UMa 
stars in the direction of one FRB, and a periodically vari¬ 
able bright blue star toward another FRB. We also revisit 
the statistical analysis of FRB110220 and FRB140514 and 
conclude that they are from one and the same repeating 
source, rather than fr om two unrelated sources as claimed by 
IPetroff et al.1 (l2015al l. The factor-of-2 decrease over 3 years 
in DM between the two bursts cannot arise in the intergalac¬ 
tic medium, invalidating any extragalactic model for FRBs 
in which the DM is a measure of the distance to the source. 
Finally, we re-assess some arguments that have been put 
forward against the Galactic flare-star scenario, and show 
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The case for nearby FRBs 3 



Figure 1. Phase-folded V and I light curves from ASAS for the 
variable star, ASAS220101-4447.1, a possible W UMa-type source 
for FRB110627. 



Figure 2. Phase-folded V and I light curves from ASAS for 
the bright variable blue star, HD174628, within the beam of 
FRB010621. 


that they are incorrect from a theoretical perspective, from 
a demonstrative observational point of view, or both. 


2 OPTICAL OBSERVATIONS AND ANALYSIS 
OF FRB FIELDS 

Following the methodology of iLoeb et al. I (l2014h to search 
for bright variable stars within the radio beams of past 
FRBs, we have analyzed data at tho se locations from 
the All-Sky Automated Survey (ASAS) lIPoimanskilllQQTll . 
ASAS uses a set of small telecopes in Chile to monitor about 
3/4 of the sky to 14 mag in the optical {V and I), with data 
for some fields going back two decades. The ASAS catalog 
also permits us to find the sky density of specific types of 
variable stars at the Galactic longitude and latitude, {I, b) 
of each FRB, and thus to estimate the chance probability of 
finding such a star within an FRB beam. 


2.1 FRB110703 

In the field of this FRB, found bv iThornton et ^ ll2013ll . 
iLoeb et al. I (120141 ') discovered, based on imaging with the 
Wise Observatory Im telescope, a V = 13.6 mag, V — I = 
0.85 mag, W UMa system with a period of P = 7.84 hr 
within the 14 arcmin diameter Parkes beam, 4 arcmin from 
the beam center. From the number (one) of W UMa’s that 
are as bright which they found outside the beam region, 
and also from th e known volume density of W UMa’s, 
iLoeb et al. I (l2014l l roughly estimated a ~ 1.7% chance prob¬ 
ability for finding such a system within the FRB beam re¬ 
gion. ASAS data for this field confirm the V and I light- 
curve shapes, amplitudes, and period for this W UMa sys¬ 
tem (named ASAS233003-0248.3 in the ASAS database). 
Furthermore, from the ASAS catalog, we hnd that the ac¬ 
tual sky density of W UMa’s in the general direction of this 
FRB, to U = 14 mag, is 0.11 deg“^, implying only a 0.5% 
random probability to find a W UMa within the 0.04 deg~^ 
Parkes beam (i.e. lower than the rough estimate, above, by 
iLoeb et al. I (l2014ll .~) 


2.2 FRB110627 

Near the position of FRBl 10627 dThornton et al.ll2013ll . we 
have found another bright {V = 12.5 mag) W UMa sys¬ 
tem, ASAS220101-4447.1. Figure [1] shows the V and I light 
curves, folded with a period of P = 10.75 hr. The period, 
the color (U — / = 0.7), and the light-curve shape (with 
amplitude 0.4 mag) are all characteristic of a W UMa star. 
The star is 29 arcmin west, i.e. about 2 beam diameters, 
from th e center of beam No. 12 in which the FRB was 
detecte dKeane fc Petrofj ll2015ll . However, beam No. 12 is 
on the outer ring of beams of the Parkes multi-beam re- 
ceiveiQ. From the technical data for FRBs at the Reasearch 
Data Australia PortaQ, the receiver orientation on the sky 
at the time of the observation was such that beam No. 12 
was almost exactly on the western side of the outer ring 
(i.e. at the “3 o’clock” position). As such, the detection in 
this single beam could have been triggered by a side lobe 
of a bright FRB at the west-offset position of the W UMa. 
The FRB/W-UMa association is thus possible in this case. 
From ASAS, the W UMa sky density in this FRB’s direction 
is 0.16 deg“^. Considering that the area within which this 
type of side-lobe detection in beam No. 12 could have been 
triggered (an annular sector of opening angle 60°, centred 
on the receiver centre, with inner radius 65 arcmin and outer 
radius 80 arcmin, i.e. at the location of the W UMa star) is 
about 0.3 deg^, the chance probability for the presence of a 
W UMa there is about 5%. 


2.3 FRB010621 

In the field o f this FRB discovered bv (iKeane et al.ll2012l l. 
ASAS data (jPoimanski fc Macieiewskil l2005| . Fig. [2] see) 
show that the bright blue star {V = 9.1 mag, V — I = 0.1), 
HD174628, 6.7 arcmin from the reported FRB radio beam 
center, is variable with a period of 2.5 days and a U-band 
amplitude of ~ 0.07 mag. Figure [3] shows the star’s spec¬ 
trum (see below) with its best-fitting model template, a 


^ see e.g., http://www.atnf.csiro.au/research/multibeam/instrument/descripti 
^ http://http://supercomputing.swin.edu.au/data-sharing-cluster/parkes-fr 
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Figure 3. Continuum-normalized Wise eShel spectrum of the 
periodic variable HD174628, with its best-fitting model stellar 
template, a B9 III star. Ha and H/3 are the prominent lines. 


B9 III giant with effective temperature T ~ lljOGOif, sur¬ 
face gravity log g = 3.5±0.5, and projected rotational broad¬ 
ening Kot sini = 70± 10 km s“^. The spectral type, period, 
and rotational broadening are all characteristic of a class of 
non-radiall y pulsating st ars called slowly pulsating B-stars 
(SPBs; e.g. iDe Catll2007l b HD174628 is somewhat unusual, 
however, in that SPBs typically have pulse amplitudes of 
only 0.01, and 0.025 at most. 

To test for radial-velocity variations, or for the possi¬ 
bility that the photometric variations are caused by some 
other star along the line of sight, whether or not physi¬ 
cally associated with the B star, we obtained high-resolution 
spectra at 8 separate epochs, as follows. Spectra with res¬ 
olution R « 10, 000 were obtained using the Wise Obser¬ 
vatory Im telescope with the eShel high resolution fiber- 
fed echelle spectrograph on four nights over a 32-day pe¬ 
riod, as detailed in Table [T] A one-hour exposure was taken 
each night. A month later, the L evy echelle spect rometer 
on the Automated Planet Finder JVogt et al.ll20lj . APF;) 
2.4m telescope at Lick Observatory was used to further ob¬ 
tain R « 100, 000 spectra at three epochs over 100 days. 
Each epoch consisted of four consecutive 15-min exposures 
that were analyzed individually. Both data sets were reduced 
and calibrated using standard procedures. 

We have processed each individual spectrum through 
the TODCOR two- dimensional correlation program 
dZucker fc Mazehlll994l l. attempting to fit each epoch with 
two stellar templates with variable radial velocities (RVs). 
In no cases do we find evidence of a secondary, late-type 
stellar component in the spectra. The spectral features in 
the spectrum consist of only a handful of broad absorption 
lines, and therefore the RV precision is not much better 


than c/R. In the Wise data, we estimate a single-epoch 
RV precision of ~ 4 km s“^. In the Lick-APF spectra, 
from the scatter in RV among each set of four individual 
exposures taken over an hour, we estimate the RV precision 
at 1 km s“^. As listed in Table [1] on all but one of the 8 
nights that the star was observed, the RV is consistent with 
a constant barycenter-corrected value of « —16 km s“^. 
One APF epoch shows a formally significant deviation to 
an RV of « —22 km s“^, but given the consistent velocities 
at all other epochs, at various phases of the 2.5 day 
photometric period, we suspect this single deviant point is 
a fluke. We detect also no significant differences in Kot sini 
among epochs. 

We have thus found no explanation for the 7%-level, 
few-day-period photometric modulation of this star, other 
than it being an SPB star with an unusually high ampli¬ 
tude. A variable line-of sight late-type star responsible for 
the variations would have contributed at least 10% of the B- 
star flux, which would have been detected in the spectrum 
and can thus be ruled out. A contact binary composed, e.g., 
of two IOM0 B-type giants with a separation of 257?© would 
have an orbital period of 3.3 days. If viewed nearly pole-on 
to the orbital plane, the photometric variability at the ±5% 
level could be due to ellipsoidal tidal distortion of the com¬ 
ponents. However, to satisfy the constraint of RV variations 
with amplitude K < 2 km s“^, despite the 200 km s“^ or¬ 
bital velocity, the system would need to be within 0.5 degree 
of pole-on. Furthermore, one would then expect the synchro¬ 
nized projected rotation velocity to be equally small, rather 
than the observed usini = 70 ± 10 km s“^. The same prob¬ 
lems would arise, but would be less severe, if the binary 
companion in the low-inclination orbit were a roughly solar- 
mass object—a main-sequence star, white dwarf, neutron 
star or black hole. For periods of 3.3 days, the B-star’s or¬ 
bital velocity would be of order 30-40 km s“^, and thus an 
inclination within 5 degrees of pole-on would be required, in 
order to satisfy the observed ARV< 4 km s“^ limits. About 
one-half of SPB stars are magnetic, with fields of order hun¬ 
dreds of Gauss and sometimes close to ~ 1 kG (iHubrig et al.l 
l2009ll . not unlike the ~ 1 kG fields of low-mass flare stars. 
B stars are relatively rare, and more so SPB stars. For ex¬ 
ample, in the 105 deg^ Kepler field, there are only about 10 
B stars brighter than 9 mag, and only 4 of them are SPBs 
(all with amplitudes 1-3 orders of magnit ude smaller than 
HD174628). A density of < 0.04 deg“^ (iMcNamara et al] 
l2012h . can thus serve as a firm upper limit for SPBs with the 
large HD174628 pulsation amplitude. The presence of this 
unusual star in the 0.042 deg^ Parkes beam in one out of 8 
FRB fields studied (random probability < 1.4%) is therefore 
intriguing. However, again, care must be taken with such a 
posteriori statistical conclusions, and therefore we cannot 
claim a secure association between the star and the FRB. 

W e note again that, for this FRB. [Bannister fc Madseiil 
(I2OI4II have used velocity-resolved observations of Ha and 
H/1 Balmer line emission from diffuse ionized gas in the FRB 
direction, combined with a Galactic rotation model, to ob¬ 
tain the emission measure of the Milky Way ionized gas 
column along the line of sight. From the emission measure, 
they estimated the Galactic electron column density. They 
concluded that the DM of the FRB, if all due to the Galac¬ 
tic interstellar medium, places the FRB source within our 
Galaxy, at a distance of 8-20 kpc. The B star HD174628, 
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discussed above, is at a distance of 0.4 — 1 kpc, depending 
on its exact spectral class and the line-of-siglit extinction 
that it undergoes. 


2.4 Non-detections 


Apart from the above three FRB fields, we have searched 
another five FRB fields for bright short-period (< 1 day) 
variable stars, but have not found any, as follows. The fields 
of FRB110220 and FRB120127 were already monitored and 
reported bv iLoeb et al. I (120141 '). who found no short-period 
stars brighter than 16 mag. This null result is confirmed 
by the shallower but longer-ter m ASAS data for these two 
fields , as well as for FRB011025 llBurke-Spohaor_&^^Barmisterl 
[20 i 3). Near the position of FRB131104 ( Ravi et al.ll20ljT 
ASAS data reveal a W UMa with V = 13.5, V — I = 1, and 
a period of P = 7.6 hr. However, the star is 28 arcmin south 
of the center of beam No. 5 in which the FRB was detected . 
Given the receiver orientation at the time (iRavi et al 
an FRB from the W UMa position would have triggered 
beam No. 12, rather than No. 5. This W UMa star is there¬ 
fore not the source of this FRB. _ 

FRB010824, the first (lUorimer et al.l 120071 ') FRB dis¬ 
covered, adjacent to the Small Magellanic Cloud (SMC), 
and at 30 Jy still the brightest one known, is also the 
only one detected in more than one focal-plane beam of 
the Parkes multi-beam receiver. The three beams with de- 
tect ions permit a more precise l ocalization of the FRB 
(^see lKulkarni et ^ (120141 '). and also lKeane fc Petrofll (l2015ll . 
who report a weak detection even in a fourth beam), to 
a roughly rectangular region about 9 arcmin by 1 arcmin 
in size. ASAS data for this area again reveal no candi¬ 
date short-period optical variables to 14 mag. This SMC 
field has also been monitored for four years by the Opti¬ 
cal Gravitational Tensing Experiment (OGLE -TV) using the 
1.3 m Warsaw University telescope in Chile (lUdalski et al.l 
I 2 OI 5 I T Within the FRB error region, we find no short-period 
variable stars in OGLE-I’V to V = 18.5 mag. From the 
OGLE-I'V data in the SMC region (iPawlak et ahll^lSl T we 
have counted the number of periodic variables redder than 
V — I = 0.5 with periods less than 1 day, to U = 18.5 mag, 
and find a density of 0.5 deg“^. These variables should in¬ 
clude all main-sequence contact binaries, including W UMa 
types. Despite the color cut, this estimate might be contam¬ 
inated by reddened main-sequence B-type contact binaries 
in the SMC itself, and therefore the density is a conserva¬ 
tive upper limit. The product of the density and the error 
area gives a chance probability of 0.13% of finding a W UMa 

within the error region. _ 

The field of FRB121102, discovered by ISnitler et al.l 
(l2014ll with Arecibo, is north of the -1-28° declination limit 
of ASAS, and in the mid-plane of the Galaxy (latitude 
b — —0.2°), where crowding and extinction make the type 
of optical followup decribed here difficult. 


2.5 Association probability of FRBs and 
W UMa’s 

Table[2]summarises some of the parameters of the ten known 
FRBs and the candidate stellar counterparts that we have 
found for some of them. In view of the above results, we re¬ 
visit the question of the probability that FRBs and W UMa’s 


are associated. In 8 FRB fields that we have searched, we 
have found one or two possible associations with a W UMa 
system, as follows. Within the be am of FRBl 1 0703. there is 
the W UMa system discovered bv iLoeb et al. I (l2014h . with a 
0.5% random probability for such a system within the FRB 
beam. However, as already mentioned, to avoid a posteri¬ 
ori statistics, this field and this system perhaps should not 
be counted. In FRBl 10627, there is a W UMa separated by 
29 arcmin from the FRB beam center, but a side-lobe detec¬ 
tion of a source at that position is a real possibility, given 
the receiver orientation at the time, with a 5% estimated 
chance probability for a W UMa. Apart from these two W 
UMa stars, in the field of FRB010621 we have found a likely 
SPB pulsator. We have not been able to show that it is con¬ 
nected to variability in a low-mass star, and it is unclear 
whether or not the SPB star itself is related to the FRB. 
We will therefore not consider it further in the accounting 
below. 

As seen in Table 1, five of the eight FRB fields have an 
a priori random probability of 0.5 to 0.8% of hosting a W 
UMa system, two fields have a probability of 4-5%, and one 
field (the SMC field, which is well localized) has a low prob¬ 
ability of 0.13%. To estimate the random probability of the 
experimental result, we find the binomial probabilities of the 
observed number of successes, or more, in each of the three 
sets of fields with similar probabilities (for the 0.5-0.8% fields 
we conservatively assume 0.8%). We then multiply the three 
binomial probabilities. The associations between FRBs and 
variable stars may exist in: zero out of 7 cases (if we do not 
count the first case, of FRBl 10703, and none of the other 
associations is real); in 1 of 7 cases (if the FRB 110628 asso¬ 
ciation is correct, and then the random binomial probability 
is 9%); or 2 of 8 cases (if we do count the case of FRB110703, 
in which case the random probability of the experimental re¬ 
sult is 3.6%.) Finally, if our accounting includes all types of 
bright short-period variable and potentially flaring stars, so 
that we include also the SPB star HD174628 in the field of 
FRB010621, then the chance probability for the experimen¬ 
tal result is (1.33 — 5) x 10“®. Given this range, we cannot 
yet claim to have established an association between the 
two phenomena. In any case, the fact that bright W UMa 
stars, or bright flaring stars of any type, were not found 
in most of the fields, does not rule out the flare-star/FRB 
connection hypothesis. Stellar sources of FRBs will likely 
have a distribution of optical luminosities and fluxes, possi¬ 
bly uncorrelated with the FRB luminosities, and therefore 
any flux-limited survey like ours for the optical counterparts 
will uncover only a fraction of the associations. 


3 FRB110220 AND FRB140514 ARE FROM 

THE SAME REPEATING SOURCE 

iPetroff et al.l (l2015al ') re-monitored the sites of a number of 
FRBs with the Parkes Radio Telescope for a total of 85 hours 
on sky, leading to the real-time discovery of FRB140514 
at the position of the previous FRB110220. Although only 
about 21 hours out of the 85 hours were spent on the field of 
FRBl 10220, all 85 hours were spent on some previous FRB 
location, and could have led to the discovery of a new FRB 
at the location of a previous one. Therefore for the pur¬ 
pose of calculating the chance probability of detecting an 
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Table 1. HD174628 spectral observations and radial velocities 


UT 

Date 

UT 

Time 

JD 

-2456700 

Observ. RV 

[km s“^] 

2014.05.14 

01:12 

91.55 

Wise —14.4 ±3.4 

2014.06.10 

23:02 

119.46 

-17.5 ±2.7 

2014.06.12 

00:58 

120.54 

-14.0 ±4.0 

2014.06.14 

00:43 

122.53 

-8.3 ±4.3 

2014.06.14 

21:07 

123.38 

-12.5 ±3.9 

2014.07.10 

06:43 

148.78 

Lick -21.5 ±1.1 

2014.10.14 

02:53 

244.62 

-15.3 ± 1.0 

2014.10.18 

03:07 

248.63 

-16.8 ±0.9 


Table 2. FRBs, candidate counterparts, and random probabilities for W UMa 


Date 

(1) 

1 

(2) 

FRB 

b 

(3) 

Ref. 

(4) 



RA(J2000) DEC 
(5) (6) 

Counterpart 
RA(J2000) DEC 
(7) (8) 

E 

deg-2 

(9) 

Prob 

% 

(10) 

Comments 

(11) 

010621 

25 

-04 

Keane et al.l [12012^ 

18 52 

-08 29 

18 51 40 

-08 24.4 

0.16 

4.4 

B giant,V = 9.1 mag, beam 10 

010824 

301 

-42 

Lorimer et al.l (12007^ 

01 18 06 

-75 12 19 



0.50* 

0.13 


011025 

357 

-20 

Burke-Snolaor & Bannisteil9 07 
12014) 

-40 37 



0.19 

0.8 


110220 

51 

-55 

Thornton et al. 

1 

201,- 

) 22 34 

-12 24 



0.10 

0.5 


110627 

356 

-42 

Thornton et al. 

1 

2013 

) 21 03 44 

-44 44 19 

21 01 01 

-44 47.0 

0.16 

5.0 

W UMa, V = 12.4 mag, beam 12 

110703 

81 

-59 

Thornton et al. 

1 

20i; 

) 23 30 

-02 52 

23 30 03 

-02 48.3 

0.10 

0.5 

W UMa, V = 13.5 mag 

120127 

49 

-66 

Thornton et al. 

1 

2013 

) 23 15 

-18 25 



0.11 

0.5 

121102 

175 

-0 

Soitler et al.l (l2 

OH) 

05 32 10 

33 05 13 

not 

searched 

0.15 

- 

Arecibo FRB, Galactic plane 

131104 

261 

-22 

Ravi et al.l ll20] 

5 

) 

06 44 10 

-51 16 40 



0.18 

0.8 

140514 

51 

-55 

Petroff et al.l ll2015a^ 

22 34 06 

-12 18 47 



- 

- 

Site of FRB110220 


Notes.- (2)(3)—Galactic longitude and latitude; (7)(8) coordinates of candidate stellar counterpart; (9) sky density, to = 14 mag, of W 
UMa stars; (10) Probability of finding a W UMa within the FRB beam area. * density of W UMa’s to U = 18.5 mag, in the SMC field. 


unrelated FRB within a re-monitored FRB region, we must 
consider all 85 hours of on-sky time. IPetroff et al.l (l2015al ) 
estimated this chance probabilty at a high level of 32%, and 
concluded that the two bursts are from different sources. 
However, based on the available information and our own 
simple and straightforward calculation, we contest this con¬ 
clusion, as detailed below. 

According to IPetroff et al.l (l2015al ). the field of 
FRB110220 was observed with a c ycle of five telescop e point¬ 
ings, called a “gridding” pattern ( Morris et al.ll20(^ . In the 
first pointing the central, 14-arcmin-diameter, beam of the 
receiver was pointed at the same position as the center of the 
beam in which FRB 110220 was detected. In the subsequent 
four pointings, the telescope was pointed 9 arcmin north, 
south, east, and west of the first pointing, and then the 
cycle was repeated. The real-time FRB140514 was detected 
only in the central beam of one of those 9-arcmin north-offset 
pointings. (Actually, from comparison of the published coor¬ 
dinates, it appears the offset was o nly about 5 arcm i n nor th 
of the FRB110220 coordinates in iThornton et ^ (l2013l) .) 
At any given moment in the real-time survey, the central 
beam was the only one covering part or all of the beam area 
of the original burst. Thus, for the purpose of estimating the 
chance probability of an unrelated burst occurring as close to 
the location of of the original burst as was, in fact, observed 
for the real-time burst, the monitored solid angle that needs 
to be considered is just the solid angle covered by the central 
beam. This solid angle is 7r(7/60)^ deg^ = 0.042 deg^. For 


an as summed all-sky FRB rate of 33 00 dav~^ dRane et al.l 
l2015l) . and the 85 hr on-sky time of the lPetroff et al.r ( 2015ah 
survey, one then expects 0.0125 of an event within the cen¬ 
tral beam. However proper account must also be taken of the 
fact that, during the entire survey, an FRB was not detected 
in any of the other 12 beams of the receiver, within which 
we would have expected 0.15 of an event outside the central 
beam, with a complementary probability for non-detection 
of 0.85. The probability for an unrelated FRB being de¬ 
tected in the central beam of the real-time experiment, and 
only there, is therefore 0.0125 x 0.85 = 1%. Conversely, we 
can state with 99% confidence that the original FRB and 
the real-time FRB are from the same, repeating, source. We 
note that, if we use the somewhat lo w er FRB rate estimates, 
2500 day“^ iKeane fc Petrofj (|2015l); Macquart fc J ohnstonI 
(l2015l) or 2000 dav~^( Burke-Spolaor fc Bannisteill2014h . the 
chance probability will go down correspondingly and the 
confidence of our conclusion rises further. 

The DMs of the original and the real-time FRBs were 
945 and 562 pc cm“^, respectively, i.e. a difference of almost 
a factor of 2. If, as we have argued, the two FRBs are from 
the same source, such a change in the intervening intergalac- 
tic free electron column density over the course of about 3 
years is impossible. This then rules out intergalactic disper¬ 
sion as the origin of the large DMs of FRBs. Conversely, 
in the flare-star scenario, there is no particular reason why 
the plasma blanket around a flare star would not change be¬ 
tween bursts separated by years. Such changes are even to 
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The case for nearby FRBs 7 


be expected, given t hat the coronal pla sma should be out¬ 
flowing, as argued bv iLoeb et al. I (l2014l i. A repeating FRB 
in the same beam with changing DM could still arise in an 
extragalactic scenario, but only if the DM is produ ced by 
electrons local to the source (e.g. I Connor et alJ (l2015h 'l. The 
source could be a true repeating one (i.e. the source is not 
destroyed in a burst), or there could be several separate FRB 
sources in an “FRB-active” galaxy, or cluster of galaxies. 


4 THEORETICAL CONSIDERATIONS 

REGARDING THE FLARE-STAR MODEL 


We now address a number of arguments that have 
been raised against the flare-star hypothesis of FRBs. 
One claim has been that the coronal plasma blanket 
around a flaring star would necessarily be too dense 
to produce the observed strict f~^ fr equency depen¬ 
dence of the p ulse arrival times of FRBs (iDenniso nl l20l4 
iTuntsovI I 2 OI 4 I 1 . For the f~^ behavior, a density rie < 
10® cm“® is required since the plas ma frequency is t/p = 
(ujy/2 n) = 0.9 GHzfw,./10^°c m~®~)^^® iTR vbicki fc LightmanI 


Il98d ). However, for example. ICetman et al.l i 2008al H) have 


found that for the brightest X-ray flares from pre-main- 
sequence stars, the coronas have electron column densities 
in the range Ne = cm“® (as implied by the DMs 

of FRBs) but also length scales as high as i? cm, 

indicating mean densities of rie ~ 10®“® cm“®. As such, one 
would expect to see the /“® frequency s weep in stellar ra¬ 
dio fla res. Indeed, there are such examples. lOsten fc BastiarJ 
(l2008ll show high-time-resolution radio data for the active 
young M-star AD Leonis, with a clear delay in signal arrival 
time with decreasing frequency. The typical drift rates are 2 
GHz s“^. Although, in this specific star, the observed drift 
rate is an order of magnitude larger than seen in FRBs, and 
it is hard to say what is the exact frequency dependence of 
the sweep, it shows that known stellar flares do sweep in 
radio frequency, at least roughly in the required sense. We 
encourage radio observers of stellar flares to obtain more 
such high-time-resolution data of more stars, to broaden the 
known phenomenology of these effects. 

Another argument against a stellar-flare origin for FRBs 
has been that radio emission would be suppressed through 
free-free absorption in the corona that produced its DM if 
the corona is limited by the radius of the star, R^,, and 
its temperatu re is limited by the viria l temperature of the 
star, ~ lO^K llLuan fc Goldreidil 120141 '). However, the opti¬ 
cal depth for free-free absorption at a frequency of 1.4 GHz 
declines with incre asing temperature T and size R of the 
emitting region as jKulkarni et al.l[2014h . 


Tf[ = 0.1 


/ ^ \ — 3/2 

VIO® K ) 


f R 


Viol® 


/ DM y 

V 10 ® cm“® pc ) 


( 1 ) 

and becomes negligible for the observed D Ms of FRBs at 
T > 10® K and R 3> lOR*. Observationally, iGetman et al.l 
((20083) I Get man et al.l (l2008bl ) show that pre-main- 
sequence flare stars often have flares with temperatures 
T ~ 10®“®K over a spatial scale R ~ 10®R* ~ lOi® cm, in¬ 
validating the argument for significant free-free absorption. 
The free-free optical depth scales with frequency as and 
therefore for the above parameters of real stellar flares, we 


would not expect to detect FRBs at frequencies much lower 
than 1.4 GHz. This is consistent with a r ecent strong limit on 
FRB d etection at 145 MHz obtained bv iKarastergiou et al. I 
ll2015h using the LOFAR array, and requiring that FRBs 
have a rising flux density spectrum between these two fre- 
que ncies. _ _ 

iKulkarni et al.l ll2014h have argued that radio emission 
from flare stars would necessarily be accompanied by bright 
X-ray flares. Considering the rate of FRBs and comparing 
it to the known statistics of all-sky X-ray variability, they 
concluded that there is a shortage of known X-ray flares, 
and hence FRBs cannot originate from flaring stars. Flar¬ 
ing stars indeed erupt in both radio and X-ray wavelengths. 
There is a well-known correlation, over many orders of mag¬ 
nitude in luminosity, between the radio and X-ray luminosi¬ 
ties of flares from stars with active coronae, whether be¬ 
tween th e peak luminositi es or the time-averaged luminosi¬ 
ties ie.g. iBenz et al.lll994h . However, there is no one-to-one 
correspondence between individual radio and X-ray flares, 
as would be required in order to derive an expected rate 
of X-ray flares from the observed FRB rate. For example, 
IWilliams et al.l (l2015l l have monitored a close M-star binary 
simultaneously in many bands, including radio and X-rays, 
and report a general lack of correlation between flares in any 
two bands. 

Thus, none of the above arguments against the flare-star 
origin of FRBs is valid, mainly because the observed prop¬ 
erties of some stellar flares already show that radio flares 
with characteristics similar to FRBs are possible. A recent 
addition to this list of similaritie s is the high circular polar¬ 
ization measured in FRB 140514 llPetroff et al.|[^15al l: radio 
bursts from flaring st ars ar e often 100% circularly polarized. 

Most recently, iKa t3 ll2015h showed that the non¬ 
monotonic dependence of burst widths on dispersion mea¬ 
sure excludes the intergalactic medium as the location of 
scattering that broadens the FRBs in time. This argues in 
favor of a local origin for FRBs, where the observed DM 
is intrinsic to the FRB sources, as expected in the case of 
flaring Galactic stars. 


5 CONCLUSIONS 

The ten published detections of FRBs have triggered a 
wealth of proposed interpretations, ranging from exotic pro¬ 
cesses at cos mological dis t ances to atmospheric and terres¬ 
trial sources. iLoeb et al. I (120141 ') have previously suggested 
that FRB sources could be nearby flare stars, and pointed 
out the presence of a W-UMa-type contact binary within 
the beam of one out of three FRB fields that they examined. 
Using time-domain optical photometry and spectroscopy, we 
now And a possible W UMa counterpart to another FRB, 
and a rare slow-pulsating B star in the beam of a third 
FRB. The random probabilities for having one or two real W 
UMa’s within the eight studied FRB beam areas are 9% 
and 4%, respectively. 

Contrary to previous claims, we conclude with 99% 
confidence that two FRBs which were discovered 3 years 
apart within the same radio beam are from the same repeat¬ 
ing source. The different DMs of the two bursts then rule 
out a cosmological origin for the DMs, but are consistent 
with the flare-star scenario with a varying plasma blanket 


© 0000 RAS, MNRAS 000, 000-000 



































































8 D. Maoz et al. 


between bursts. Finally, we have shown that the theoreti¬ 
cal objections that were raised against a Galactic origin of 
FRBs are incorrect because they are circumve nted by the ob¬ 
served properties of some stellar radio flares (iGetman et al.l 
l2008al la). We conclude that the flare-star-origin hypothesis 
for FRBs is still a strong player in the game, and is con¬ 
sistent with all current observations. In contrast, the idea 
that the intergalactic medium is responsible for the DM of 
FRB sources is in direct contradiction to the observation of 
a changing DM between two FRBs that, as we have argued, 
are almost certainly from the same repeating source. 
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